Exploiting quantum coherence of polaritons for ultra sensitive detectors 
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Besides being supernuids, microcavity exciton-polariton condensates are capable of spontaneous 
pattern formation due to their forced-dissipative dynamics. Their macroscopic and easily detectable 
response to small perturbations can be exploited to create sensitive devices. We show how controlled 
pumping in the presence of a peak-dip shaped potential can be used to detect small externally applied 
velocities which lead to the formation of traveling holes in one dimension and vortex pairs in two 
dimensions. Combining an annulus geometry with a weak link, the set up that we describe can be 
used to create a sensitive polariton gyroscope. 
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Microcavity exciton-polaritons are two-dimensional 
half-light half-matter quasi-particles that result from the 
hybridisation of quantum well excitons and photons in a 
planar Fabry- Perot resonator. At low enough densities, 
they behave as bosons and if the temperature is lower 
than some critical value they undergo Bose-Einstein con- 
densation [TJ H] • Recent experiments have investigated 
exciton-polariton condensation and the phenomena asso- 
ciated with it, such as pattern formation [3j[4], quantised 
vortices and solitons [5 , 6], increased coherence and the 
cross-over to regular lasing; recent reviews of the field 
can be found in [7, 8J. These discoveries have opened 
a path for new technological developments in optoelec- 
tronics (ultra-fast optical switches, quantum circuits), 
medicine (compact terahertz lasers) and power produc- 
tion (hybrid organic-inorganic solar cells). However one 
fundamental property of exciton-polariton condensates, 
their ability to flow without friction at moderate veloc- 
ities |9j [10] , has never been considered for technological 
use. The existence of frictionless flow in polariton con- 
densates sets them apart from other solid-state quantum 
systems and links them to cooperative fluids, such as 
superfluid helium or atomic condensates. In this letter 
we propose an idea that may lead to a new generation 
of ultra-sensitive devices based on the quantum coher- 
ence and superfluidity of exciton-polariton condensates. 
The important ingredients of this novel technology come 
from the unique properties of polariton condensates: (1) 
they are non-equilibrium systems capable of pattern for- 
mation; (2) their dynamics is controlled by the balance 
between gain, due to continuous pumping, losses and 
non-linearity; (3) polaritons condense at relatively high 
(even room) temperature thanks to their very small ef- 
fective mass; (4) one can easily engineer any external 
landscape and vary pumping in space and time; (5) po- 
lariton condensates form quantised vortices in response 
to slight changes in the environment: when flow exceeds 
the critical velocity, when fluxes interact, when pumping 
powers exceed a threshold for pattern forming instabil- 



ities, when the magnetic field exceeds a threshold, etc. 
These properties allow one to prepare the system in a 
state slightly below the criticality for vortex formation, 
so that a tiny external perturbation will take it over the 
criticality, leading to a macroscopic and easily detectable 
response. Similar principles underpin superconducting 
single photon detectors, biased just below their transi- 
tion temperature. Therefore, we propose exploiting su- 
perfluid properties and sensitivity to vortex formation of 
polariton condensates to create sensitive devices that re- 
spond to slight changes in rotation rate, gravity or mag- 
netic field. We will illustrate the main principles using 
the example of an exciton-polariton gyroscope. 

Superfluid helium gyroscopes have already been shown 
to provide sensitive means for detecting absolute rota- 
tions [TTJ [12] , but required very low temperature for op- 
eration. A superfluid in an annulus (ring) that is slowly 




FIG. 1: (color online) Schematic of the experiment: pumped 
polariton condensate channel with a weak link. If the flow 
across the link is high enough, vortices are generated. 

rotated remains motionless as there is no friction with 
boundaries. But if a partition with a small opening (weak 
link) is inserted in the annulus, a flow in the direction op- 
posite to the rotation will be generated through it. To 
leading order [12 , the velocity U across the opening is 



2 



related to the angular speed of rotation Q by the relation 
U = VtR 2 /5, where R (6) is the width of the annulus 
(opening). When U exceeds some critical value, which is 
of the order of the Landau critical velocity U c , vortices 
are generated and detected as phase slips. For the case of 
superfTuid He4, U c ~ 50m/s,R ~ 0.1m, 5 ~ 100 nm, al- 
lowing accurate detection of the Earth's daily rotation 
rate. It was also established that superfluid rotation 
sensors, similar to atomic beam gyroscopes, belong to 
the same class of quantum interference effects as Sagnac 
light- wave experiments [13]. The same principles can be 
applied to a polariton gyroscope, schematic of which is 
given in Fig. [TJ which offers the advantage of potentially 
operating at room temperatures. The critical velocity for 
vortex formation in polariton condensates [10 is three- 
four orders of magnitude higher than that in superfluid 
He4, a fact that would reduce the sensitivity of the po- 
lariton gyroscope. However, the pattern forming proper- 
ties of polariton condensates allow one to severely reduce 
the velocities needed for the vortex formation to occur 
and to use pumping intensity as the control parameter. 
We propose to insert a peak-dip shaped potential, which 
will further accelerate the flow, at the position of the 
weak link. Such a potential can be prepared either us- 
ing a combination of etching [TIHTS] and stress induced 
traps [1 , or directly defining blue-shifted trap potentials 
via spatial light modulators [3]. Even in the absence of 
any externally applied flow, velocity fluxes connecting re- 
gions where the density is low (at the potential peak) to 
regions where the density is high (at the potential dip) 
will be generated; in between them there will be a point 
where the condensate speed takes its maximum value. 
Close to such a point, the density has a local minimum 
(Bernoulli effect). A slight change of external conditions, 
e.g. an increase of the potential strength, or a decrease 
of the pumping intensity, will further lower the density. 
If the perturbation is large enough, the density minimum 
will reach zero and a vortex pair will be emitted, as we 
illustrate below using a mean-field model of polariton 
condensates. The scheme for measuring rotations then 
becomes quite straightforward: having prepared a po- 
tential such that the system is in a slightly subcritical 
configuration for the range of rotation speeds that one 
intends to measure, the strength of the pumping inten- 
sity is varied, recording the moment when vortices start 
to nucleate. From this one deduces the back- flow through 
the aperture and, therefore, the rotational velocity. 

Polariton condensates at temperatures much lower 
than the critical temperature for condensation can be ef- 
fectively described using the complex Ginzburg-Landau 
equation (cGL) [I9j[20]. If a condensate is flowing with 
bulk speed U along the x axis, the cGL in a frame co- 
moving with the condensate takes the form 
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where V(x,y) describes an external potential, a is a 
pump rate, a is a nonlinearity which causes pumping 
to reduce as density increases and r] ~ 0.1 is an energy 
relaxation term [21] which causes the evolution towards 
states of lower energy. Linear losses are included in a. 
Eq. ([I]) is stated in harmonic oscillator units, assuming 
that V = uor 2 /2 and measuring energy in units of Hlj/2, 
length in units of the oscillator length I = ^h/muj, and 
time in units of where uo is the oscillator frequency 
of the trapping potential. For m and Huj we take the 
values found by Balili et al [TJ [23]: Hlj = 0.066 meV, 
m = 7 x 10 _5 m e and for a and a we follow estimates 
given in Refs. [19j[22]: ° = 0.3, < a < 10. 

In order to illustrate the main principles of the vor- 
tex formation mechanism, we consider first the simpler 
ID case, with an external potential of the form V(x) = 
Vo (exp (— (x — xo) 2 ) — exp(— (x + xo) 2 )). The peak-dip 
shape of this potential allows for a better acceleration of 
the condensate than a single bump; moreover, in 2D the 
dip has also the useful role of temporarily trapping vor- 
tices and making their detection easier. After perform- 
ing a Madelung transformation setting ip = v /pexp(z^>), 
equation in a frame at rest (U = 0) and for a station- 
ary state becomes the system 



d x (pu) = (a- ap-r]p)p 

p = p + u 2 + V(x)-d 2 xy fp/,fp 
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where p is the chemical potential. The speed u = V0 
takes its maximum at some point between — xq and xo, 
the maximum and minimum of the potential. It follows 
from Eq. (|3|, which is a generalized version of Bernoulli's 
equation, that close to the maximum of the speed, den- 
sity has a local minimum. The numerical results 
displayed in Fig. [5] show that max('u) is an increasing 
function of Vq and min(p) a decreasing one. This can 
be verified analytically deriving an approximate solution 
as follows. For very small Vo it is natural to expect a 
dip-peak density profile; however, for higher values of Vb, 
there will be an additional minimum due to the Bernoulli 
effect. Therefore we take a density ansatz of the form 
p =a/(a + rj) + B exp(—b(x — x 3 ) 2 ) + Cexp(-c(x — xi) 2 ) 
—Qex.p(—q(x — X2) 2 ), where 6, B, c, C, q, Q, x±, £2, £3 are 
free parameters. Substituting in ([2| and integrating one 
has: 
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The potential V(x) is then determined from (|3| and con- 
tains the unknown parameters B, C, Q : b, c, g, xi, £2, x%. 
The values of these parameters can be fixed requiring 
V(x) to fit the original potential V(x). A comparison 
of the analytical and numerical solutions for u is shown 
in panel (a) of Fig. [2j while numerical solutions for p are 
shown in panel (b). The behavior of max(u) and min(p) 




FIG. 2: (color online) (a) Numerical (black dashed line) and 
analytical estimates (red solid line) of velocity profiles as Vb 
is increased from 1 to 4. (b) Numerical solutions for p with 
Vo = 1.4,2.0,2.6,3.2,3.8,4.4. Higher values of Vo correspond 
to lower values of min(p). In both figures xq = 1.5. 
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FIG. 3: (color online) (a) ID. Profiles of max(«) (numerical: 
black dashed; analytical: red continuous line) and of min(p) 
(numerical: blue dot-dashed; analytical: red continuous line) 
as a function of Vb, with xo = 1.5. The theoretical solu- 
tion is only valid up to Vb ~ 4.0. (b) 2D. Numerical profiles 
of min(yo(x,0)) (blue dot-dashed line) and of m&x(u x (x,0)) 
(black dashed line) as a function of Vo, xo = 1.5. 

with Vo is mostly linear up to Vo = 4, see Fig. [3] panel 
(a). A deviation from the linear regime leads quickly, for 
Vo > 4.4, to the loss of stability of stationary solutions. 
The time dependent solutions found for Vb > 4.4 are 
characterized by the periodic emission of traveling holes 
[24] . which travel for a short time before dissipating, see 
Fig. [4j The density hits zero when the hole is emitted. 

In two dimensions we take a similar potential: 
V = V (exp(-((x + x ) 2 +y 2 ))~ exp(-((x - x ) 2 +y 2 )))- 
The dynamics is similar to that of the ID case: min(p) is 
again a decreasing function of Vb, with a linear behavior 
for small values of Vb, see Fig. |3]panel (b). When Vo 
exceeds a critical value, stationary solutions become 
unstable and vortex pairs form. If the system is close 
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FIG. 4: (color online) Traveling holes, identified by the local- 
ized low- density blue region, form around x = 0, and travel 
up to x ~ 1.5. The process repeats periodically in time. 
(Vo = 4.8, xo = 1.5, a = 4, U = 0). 

enough to criticality, the drop of min(p) to zero and 
consequent onset of vortex nucleation can also be started 
by a small increase of the external velocity U or decrease 
of the pumping strength a. The process is shown in 
Fig. [5j where the onset of vortex nucleation is due to an 
increase of U. Representatives of critical and subcritical 
states are shown in Fig. [6j There is, however, one 
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FIG. 5: (color online) Real space plots of the zero level of 
real (black continuous line) and imaginary (grey dashed line) 
part of superimposed over the contour plot of p. Regions of 
higher density are red. (a) A slightly subcritical configuration 
(a = 4, Vo = 16.8, x = 1.5, U = 0); (b) A slightly overcritical 
configuration (a = 4, Vb = 16.8, xo = 1.5, U = 0.2) at the 
moment of vortex pair formation. A red arrow shows the 
direction of the flow; a red dot shows the point where the 
vortex pair is being generated. 

notable difference with the ID case: while in the latter 
the emission of traveling holes sets in at a finite value 
of min(p), in 2D there are stationary solutions all the 
way down to min(p)= 0, see the right panel of Fig. [3] 
An example calibration curve is shown in Fig. [7| where 
the critical values of a which start vortex nucleation 
are recorded for fixed Vo and different values of U. 
The boundary U(a) between subcritical and critical 
behavior is, for a bulk speed U < 0.5, nearly linear, 
hence the sensitivity of the apparatus scales linearly 
with a or, equivalently, with the bulk density a /(a + rj) 
of the condensate. In the numerical simulations that we 
performed, the system showed a clear transition between 
subcritical and critical states at velocities of the order of 
10 -4 of the critical velocity U c = \J p/m for this system. 



4 




FIG. 6: (color online), (a) Density profile above criticality 
(Vb = 14, xo = 1.5, a = 4, U = 0.4). The vortex pair has been 
temporarily trapped by the dip in the potential, (b) Density 
profile for a stationary solution below criticality (Vb = 14, 
xo = 1.5, a = 4, [7 = 0). The 3D plots show the density 
profile, with the potential superimposed in yellow. The 2D 
contour plot gives the projection of the density values on the 
xy plane, with regions of higher density colored in red. The 
direction of the flow is from the left to the right. 
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FIG. 7: The main plot shows part of the boundary between 
critical and subcritical regions of phase space with respect to 
the parameters U and a for Vb = 14, xq — 1.5. The inset has 
been obtained with a different values of the potential, Vb = 
16.8, xo = 1.5 and gives an idea of the attainable resolution 
for small values of U. 

This puts the sensitivity of polariton gyroscopes in the 
same range as that of superfluid helium gyroscopes, 
but at the advantage of potentially operating at room 
temperature. 

In summary, we proposed an idea for creating polari- 
tonic sensitive devices, such as a superfluid gyroscopes, 
based on their macroscopic response to small perturba- 
tions. We studied the mechanism of the formation of 
traveling holes in ID and vortex pairs in 2D in an exter- 
nally imposed peak-dip shaped potential when the height 
of this potential (the difference between maximum and 
minimum) exceeds the threshold. 
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